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1.  Introduction 


Currently,  there  is  great  interest  in  acoustics  within  the  Army  because  of  the  potential 
of  its  use  for  non-line-of-sight  battlefield  detection  and  location.  Turbulence  has  been 
identified  as  a  candidate  cause  for  higher  than  expected  sound  levels  within  shadow 
zones,  those  regions  where  ray  theory  would  indicate  no  sound  could  penetrate. 
Traditional  turbulence  theory  has  been  developed  around  a  model  of  an  energy  cascade 
through  a  series  of  eddies  (hereinafter  called  turbules)  of  decreasing  size  until  viscosity 
effects  dissipate  the  energy.  Traditional  turbulence  theory  is  more  completely  covered 
elsewhere  (Noble  and  Auvermann,  1993).  The  theory  goes  immediately  into  statistical 
averages  with  no  further  interest  in  the  internal  structure  of  the  turbulent  region.  This 
traditional  theory  is  therefore  termed  the  "statistical  approach"  to  turbulence.  More 
recently  as  experimental  evidence  has  accumulated,  interest  has  developed  in  the 
"structural  approach"  to  turbulence  wherein  the  effect  of  internal  structure  is  more 
adequately  accounted  for.  There  is  a  natural  parallel  between  optical  scattering  from 
aerosols  and  acoustical  scattering  from  turbules.  This  parallel  comes  from  the  similarity 
between  the  picture  of  atmospheric  aerosols  as  a  collection  of  spherical  particles  of 
varying  index-of-refraction  and  size  and  the  picture  of  a  turbulent  region  as  a  collection 
of  spherical  eddies  with  varying  wave  speeds  and  sizes.  To  determine  scattering  from 
such  collections,  the  scattering  from  each  size  of  scatterer  must  be  determined  and  then 
the  sum  from  some  specified  size  distribution  calculated.  This  is  exactly  the  function  of 
the  optical  scattering  code  AGAUS. 

Program  AGAUS  was  developed  several  years  ago  to  determine  electromagnetic 
scattering  properties  of  small  spherical  particles.  This  report  is  concerned  with  the 
conversion  of  the  AGAUS  code  from  an  optical  scattering  code  to  an  acoustical  scattering 
code.  The  notion  is  that  the  temperature  distribution  within  a  turbule  creates  an 
accompanying  sound  speed  distribution  that  will  cause  the  sound  speed  within  the  turbule 
to  be  different  from  that  outside  the  turbule.  There  is  thus  an  acoustic  index  of 
refraction  discontinuity  caused  by  the  turbule  in  the  same  way  an  aerosol  particle  causes 
an  optical  index  of  refraction  discontinuity.  Acoustical  scattering  is  also  caused  by  the 
velocity  distribution.  This  is  not  treated  here.  The  index  of  refraction  variation  within 
a  turbule  is  treated  as  an  average  quantity  so  it  is  uniform  within  the  spherical  volume 
occupied  by  the  turbule.  One  other  difference  between  optical  scattering  and  acoustic 
scattering  is  that  the  index  variations  in  turbulence  are  very  slight  as  compared  to  index 
variations  in  aerosols  distributions.  One  should  be  careful  to  calculate  with  sufficient 
precision  so  that  the  effects  are  not  lost  in  computer  round  off  errors.  Small  temperature 
variations  between  the  inside  and  outside  of  a  turbule  result  in  refractive  indices  near  1 . 
Hence,  the  conversion  of  program  AGAUS  from  single  precision  to  quadruple  precision 
was  accomplished,  resulting  in  the  code  called  QGAUS.  This  conversion  is  described 
elsewhere  (Nguyen,  1993).  The  conversion  to  acoustic  scattering  began  with  QGAUS. 

Electromagnetic  scattering  code  AGAUS  utilizes  Mie  theory  for  determining  the 
scattering  properties  of  spherical  homogeneous  particles.  Mie  theory  predicts  the 
scattering  and  absorption  of  electromagnetic  radiation  of  wavelength  X  by  spherical 
particles  of  radius  r  having  refractive  index  m  =  n  +  ik  relative  to  the  propagation 
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medium.  A  complex  refractive  index  is  used  in  the  optical  code  AGAUS  to  account  for 
a  lossy  propagation  medium,  and  this  usage  is  retained  in  the  acoustic  scattering  code 
ASCT.  Program  runs  of  ASCT  were  performed  with  the  imaginary  component  of  the 
complex  refractive  index  m  set  to  zero  as  only  real  refractive  indices  were  of  interest. 
Particles  in  space  attenuate  incident  electromagnetic  radiation  by  scattering  the  radiation 
in  different  directions  and  absorption  within  the  particle.  Extinction  is  the  sum  of  the 
scattering  and  the  absorption.  Extinction,  scattering,  and  absorption  of  a  particle  are 
characterized  by  effective  "bloc king”  areas  or  cross-sections  (Ccrt,  C^J.  The 
angular  distribution  of  the  scattered  radiation  of  spherical  particles  is  characterized  by 
the  phase  function  p(0),  a  function  of  the  single  variable  6  that  is  the  angle  between  the 
direction  of  propagation  and  the  direction  to  the  observation  position.  The  cross-sections 
and  phase  function  are  dependent  on  the  particles’  refractive  index  and  size  parameter 
x  =  2-rr/X,  where  r  is  the  particle  radius  and  X  is  the  radiation  wavelength.  Mie  theory 
provides  a  method  for  calculating  the  cross-sections  and  phase  function  as  a  series 
expansion  of  the  field  into  partial  waves  (monopole,  dipole,  quadrapole,  etc).  A  detailed 
description  of  the  theory  can  be  found  in  Bohren  (1983). 

Because  electromagnetic  and  acoustic  wave  propagation  both  obey  the  wave  equation,  the 
same  theory  can  be  used  to  determine  acoustic  scattering  properties  of  spherical  turbules. 
For  electromagnetic  waves  Mie  theory  matches  partial  wave  expansion  coefficients  at  the 
particle  boundary  so  that  both  the  electric  field  (E)  and  the  magnetic  field  (H)  are 
solutions  of  the  wave  equation  inside  and  outside  the  particle.  When  Mie  theory  is 
applied  to  acoustic  waves,  the  solution  to  the  wave  equation  involves  matching  partial 
wave  expansion  coefficients  for  the  sound  wave  pressure  and  velocity  at  the  turbule 
boundary.  A  detailed  derivation  of  expansion  coefficients  expressions  for  acoustical 
scattering  by  a  uniform  density  spherical  turbule  has  been  provided  through  a  separate 
contract  (appendix  A).  Analytical  formulae  from  this  derivation  were  used  in  the 
conversion  of  QGAUS  to  the  acoustic  scattering  code. 

Details  of  the  conversion  of  electromagnetic  scattering  code  QGAUS  to  acoustic¬ 
scattering  code  ASCT  (for  Acoustical  Scattering  from  Turbules)  are  presented  in  section 
2.  Representative  results  for  the  cross-sections  and  phase  function  are  presented  in  tables 
and  graphs  in  section  3.  Code  verification  is  obtained  by  a  comparison  of  code  results 
for  extinction  with  approximate  analytical  formula  results.  Conclusions  concerning  the 
conversion  and  results  are  presented  in  section  4.  Appendix  B  describes  the  four  size 
distributions  that  may  be  chosen  in  ASCT.  The  electromagnetic  scattering  formulas  that 
are  the  basis  of  AGAUS  are  included  in  appendix  C  to  compare  with  the  corresponding 
acoustic  formulas  more  easily. 

2.  Code  Conversion 

A  detailed  description  of  AGAUS  can  be  found  in  Miller  (1983).  Since  Program 
QGAUS  is  identical  to  Program  AGAUS  except  that  all  calculations  are  done  in 
quadruple  precision,  the  description  in  Miller  (1983)  is  directly  applicable  to  QGAUS. 
All  work  described  in  this  report  and  in  Nguyen  (1993)  was  done  on  a  Hewlett-Packard 
9000/780  mainframe  using  HP  FORTRAN  77/HP-UX.  The  only  significant 
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computational  differences  between  AGAUS  and  QGAUS  result  from  the  absence  of  quad 
precision  complex  arithmetic  operations.  Appropriate  quad  precision  real  operations 
were  substituted  where  necessary  and  combined  to  give  a  quad  precision  complex  result. 
Program  QGAUS  has  10  modes  of  operation  defined  by  10  different  particle  size 
distribution  models.  Inputs  to  QGAUS  include  wavelength,  complex  refractive  index, 
and  11  integer  program  control  parameters.  Program  output  consists  of  cross-sections 
and  phase  function  tables.  Mie  calculations  are  performed  in  subroutine  MIEGX. 
Variable  parameters  to  subroutine  MIEGX  include  the  size  parameter  x  =  2ra/X  and  the 
refractive  index  n  =  m  +  ik.  MIEGX  returns  scattering  efficiency  factors  from  which 
the  scattering  cross-sections  and  phase  function  are  calculated. 

Conversion  of  program  QGAUS  to  an  acoustic  scattering  code  involved  redefinition  of 
the  formulas  used  in  subroutine  MIEGX  for  determining  the  efficiency  factors. 
Subroutines  AGAUG  and  WATER  were  deleted  since  the  new  code  will  not  be  used  to 
treat  hygroscopic  effects.  Control  parameter  IW  was  redefined  as  ICQ  and  is  used  for 
printing  either  efficiency  factors  (ICQ  =  0)  or  cross-sections  (ICQ  =  1).  Program 
QGAUS  used  units  of  micrometers  for  wavelength  and  particle  radius  and  square 
centimeters  for  the  cross-sections.  The  modified  code  uses  units  of  meters  for  the 
wavelength  and  turbule  radius  and  square  meters  for  the  cross-sections.  Six  of  the  10 
distribution  models  were  deleted  since  they  had  no  bearing  on  the  problem  of  acoustic 
scattering  by  spherical  turbules.  A  description  of  the  four  remaining  models,  as  well  as 
the  inputs  required  by  each,  is  given  in  appendix  B. 

Subroutine  MIEGX  has  been  modified  to  compute  scattering  properties  of  spherical 
turbules  via  the  expansion  coefficient  a,,  which  corresponds  to  the  induced  electric 
multipole  moment  in  electromagnetic  scattering  problems.  The  expansion  coefficient  af 
is  defined  as  follows: 


mj[  (mxyjjx)  -  (x) 

mj',(.mx)h"Xx)  -  7,(m*)Vl,'«’ 


(1) 


where  jf(z)  and  h((z)  are  the  spherical  Bessel  function  and  the  Hankel  function  of  the 
third  kind,  respectively.  The  prime  means  "derivative  with  respect  to  the  argument." 
For  computational  purposes  a,  is  rewritten  as 


[mAfjnx)  +  -  9t45(.,«] 

[mA,(mx)  +  —]?,(*)  -  £,_,(*) 


(2) 


where 
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(3) 


Afi*  -  ^ 

'  jfinx)  jt(mx ) 


«♦!) 

ffU 


J_  +  W"1*)  _  _1_ 

mx  y((mx)  mx’ 


and  £(x)  is  the  Ricatti-Bessel  function  of  the  third  kind  and  is  related  to  the  Bessel 
function  j,(z)  by 


5,(z)  *  -  izyJLzY 


(4) 


The  difference  of  the  first  two  terms  of  the  quantity  Af(mx)  is  computed  using  Lentz’s 
continued  fraction  method.  Detailed  information  on  this  method  can  be  found  in  Lentz 
(1976). 

The  cross-sections,  scattering  amplitude,  differential  scattering  amplitude,  and  phase 
function  are  defined  as  follows: 

The  extinction  cross-section  is  given  by 


kr  to 


(5) 


the  scattering  cross-section  is  given  by 


=  iri:<2f+i>i  «.i 

kr  to 


(6) 


the  absorption  cross-section  is  given  by 


<W*)  “ 


(7) 


the  radar  cross-section  is  given  by 
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(8) 


1  £(2l+l)(-l),al|  , 

IT  t*C 


the  scattering  amplitude  is  given  by 


Am  -  4  E(2t+1*,  ?<(***)• 

**  i*« 


where  P,(cos  9)  is  the  Legendre  polynomial  function  and  is  computed  in  MIEGX  using 
the  following  recurrence  formula 

(2t  +  l)xP,(x)  =  (f  +  l)P/+,(x)  +  fPM(x), 

the  differential  scattering  cross-section  is  given  by 


o(L0)  -  |  **,6)1  2, 


and  the  phase  function  is  given  by 


m  = 


for  k  =  2t/X. 

The  radar  cross-section  and  phase  function  were  defined  according  to  those  given  in 
Bohren  (1983).  Subroutine  MIEGX  computes  efficiency  factors  Q„„  QK„  Q,bs,  and  Qr,d, 
from  which  the  cross-sections  are  computed.  The  efficiency  factors  are  multiplied  by 
the  geometrical  sphere  cross-section  to  obtain  the  true  cross-sections  C;  =  xrQ,  in  the 
main  program. 

The  efficiency  factors  are  defined  as  follows: 


rw 
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<?«.  «  4r<2<+1)!  *«t 2’ 

or  i-o 


f  13) 


4  1 

Q*  *  3 1  ECZ^iW-i)'  i  - 

JT  <-0 


(14) 


As  mentioned  above  the  difference  of  the  first  two  terms  of  the  quantity  A,(mx)  is 
computed  via  the  continued  fraction  method  (Lentz,  1976).  This  method  is  valid  for  f  = 
1,2,3,. ..etc.  Equations  (5)  through  (14)  require  computation  of  A,(mx)  for  !  -  0,  which 
is  not  possible  with  the  continued  fraction  method.  Subroutine  AO  was  developed  to 
handle  computation  of  Af(mx)  for  l  -  0.  \  is  defined  as  follows: 


v4„(mx)  * 


70(**) 


(15) 


where 


j-iimx) 


cos(mx) 
mx  ’ 


j0(w) 


sm(ntx) 

mx 


(16) 


Substituting  equation  (16)  into  equation  (15)  yields, 


AQ(mx) 


cos(mx) _ 1_ 

sin(mx)  mx 


(17) 


The  argument  mx  may  be  complex.  Hence,  computing  of  the  ratio  cos(mx)/sin(mx)  is 
done  in  subroutine  AO  by  defining  cos(mx)  and  sin(mx)  as  follows: 


cos(ru)  * 


.Imx 


sin (mx)  * 


.Imx 


2  i 


(18) 
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Once  the  quantity  A,(mx)  is  computed  for  f  =  0,  the  remaining  terms  (for  f  = 
1,2,3... etc)  are  computed  via  the  continued  fractions  method  already  coded  m  subroutine 
MIEGX. 

For  small  si2e  parameter  x,  the  dominant  term  in  the  summations  of  equations 
(5)  through  (14)  is  a,,.  Subroutine  MIEGX  bypasses  the  continued  fraction  calculations 
in  this  case,  calls  subroutine  AO,  and  computes  the  efficiency  factors  and  phase  function 
using  only  the  term  ao. 

A  comparison  of  the  formulas  implemented  in  program  AGAUS  to  those  implemented 
in  the  acoustic  code  reveals  both  similarities  and  differences  between  the  two  sets  of 
formulas.  Appendix  C  contains  the  definitions  of  the  formulas  implemented  in  program 
AGAUS.  The  corresponding  formulas  implemented  in  the  acoustic  code  are  given  by 
equations  (2),  (3),  (12),  (13),  (14),  and  (11). 

A  comparison  of  equation  (2)  with  equations  (C-l)  and  (C-2)  shows  that  the  formula  for 
the  expansion  coefficient  a,  given  by  equation  (2)  is  most  like  the  formula  for  the 
expansion  coefficient  b,  given  by  equation  (C-2).  Equation  (2)  and  (C-2)  differ  in  the 
term  (f+c)/x  where  c  =  0  in  equation  (C-2)  and  c  =  1  in  equation  (2).  Another 
difference  between  equation  (2)  and  (C-2)  is  in  the  definition  of  the  term  Af(mx). 
Equations  (3)  and  (C-3)  give,  respectively,  the  acoustic  and  electromagnetic  code 
formulas  for  the  term  A,(mx).  The  two  formulas  are  the  same  except  for  the  absence 
of  the  term  -1/mx  in  equation  (C-3).  The  formulas  for  the  extinction,  scattering,  and 
radar  efficiency  factors  for  the  acoustic  and  electromagnetic  codes  are  given  in  equations 
(12)  through  (14),  and  equations  (C4)  through  (C6),  respectively.  The  formulas  given 
by  equations  (C-4)  through  (C-6)  are  dependent  on  both  expansion  coefficients  a,  and  b, 
whereas  those  given  by  equations  (12)  through  (14)  are  dependent  only  on  the  expansion 
coefficient  af.  Also,  the  summation  coefficient  4/x2  shown  in  equations  (12)  and  (13)  is 
twice  the  coefficient  2/x2  shown  in  equations  (C-4)  and  (C-5),  and  four  times  the 
coefficient  shown  in  equation  (C-6). 

The  formulas  for  the  phase  function  are  given  by  equations  (11)  and  (C-13). 
Computation  of  the  phase  function  in  program  AGAUS  is  also  dependent  on  the 
expansion  coefficients  a,  and  b,  and  the  angular  factors  t,{9)  and  t,(6).  In  contrast,  the 
phase  function  given  by  equation  (11)  is  dependent  on  the  expansion  coefficient  a,  and 
the  Legendre  polynomial  function  P,(cos  9).  Furthermore,  the  coefficient  1  /2k2  shown 
in  equation  (C-13)  is  half  of  the  coefficient  shown  in  equation  (11).  Finally,  the 
summations  computed  in  program  AGAUS  (see  equations  (C-4)  through  (C-8)  begin  with 
t  -  1,  whereas  the  summations  computed  in  the  acoustic  code  (see  equations  (9),  and 
(12)  through  (14))  begin  with  l  -  0.  Appropriate  changes  were  made  to  the  formulas 
implemented  in  program  QGAUS  based  on  the  differences  mentioned  above  to  convert 
it  to  acoustic  scattering  code  ASCT. 

Program  ASCT  operates  under  four  distinct  modes.  In  type  0  mode  the  user  enters  a  set 
of  radii  and  corresponding  absolute  number  densities.  The  program  calculates  the  cross- 
sections  and  phase  function  for  each  radii  and  sums  the  contribution  of  each  to  get  the 
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total  cross-sections  and  phase  function.  In  type  1  mode  a  log-normal  distribution 
function  is  used  to  determine  the  radii  at  which  the  cross-sections  and  phase  function  are 
calculated.  In  type  2  mode  crusa-section  and  phase  function  are  calculated  for  a  single 
turbule  radius.  In  type  3  mode  a  table  of  cross-sections  (ICQ  =  1)  or  efficiency  factors 
(ICQ  —  0)  versus  turbule  radius  is  generated.  Appendix  B  describes  the  four  modes  as 
well  as  program  input  requirements.  The  code  verification  method  and  the  results 
obtained  are  discussed  in  section  3. 

3.  Cross-Section  and  Phase  Function  Results 

The  code  was  verified  by  comparing  the  values  for  extinction  efficiency  obtained  from 
program  ASCT  to  those  obtained  from  "anomalous  diffraction"  formula  (ADF)  given 
below, 

=  2  -  — sinp  +  -^-(1  -  cosp),  (19) 

p  p2 


where  p  =  2x  |  m-1 1 .  Equation  (19)  was  derived  in  Van  de  Hulst  (1957)  and  is  good  for 
large  size  parameter  x  (that  is,  for  X  <  r)  and  refractive  indices  near  1  (that  is.  m  =  1 
±  e).  For  small  size  parameter  x  (that  is,  X  >  r),  the  following  equation  was  derived 
for  comparison, 


<?« 


-  H 


(20) 


Equation  (20)  was  derived  following  the  method  used  in  Bohren  (1983)  for  particles 
small  compared  to  the  wavelength.  See  also  reference  Morse  (1968),  page  427. 

Program  ASCT  and  ADF  results  for  extinction  efficiency  are  shown  in  figures  1  and  2, 
for  refractive  indices  m  =  1.0000015  and  m  =  1.33333,  respectively.  The  authors 
recognized  that  1.33333  is  far  too  large  to  represent  an  acoustic  index  of  refraction. 
This  index  is  used  here  solely  as  a  means  to  exercise  and  check  the  code  in  the  large  size 
parameter  regime.  To  get  results  where  Qert  flattens  out  for  an  index  of  1.02  would 
require  an  enormous  size  parameter,  for  which  the  run  time  would  have  been  prohibitive. 
Figure  1  gives  the  extinction  curve  for  a  turbule  of  radius  0.5  m  <  r  <  105  m, 
wavelength  X  =  5  m,  and  refractive  index  m  =  1.0000015.  The  corresponding  data  for 
figure  l  is  shown  in  table  1.  Program  ASCT  extinction  efficiency  values  are  given  in 
column  Q(EXT),  and  e  extinction  predicted  by  ADF  is  given  in  the  column  so  marked. 


Figure  2  shows  the  extinction  efficiency  curve  for  a  turbule  having  refractive  index  m  = 
1.33333  and  radius  0.5  m  <  r  <  105  m.  Table  2  shows  the  corresponding  data  for 
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figure  2.  For  refractive  index  m  =  1.0000015,  the  average  percent  error  between  the 
program  and  formula  extinction  efficiencies  is  *»  0.04  percent;  and  for  refractive  index 
m  =  1.33333,  the  average  percent  error  is  »  7.4  percent. 

Figures  3  and  4  are  plots  of  the  percent  error  versus  size  parameter  x  for  refractive 
indices  m  =  1.0000015  and  m  =  1.33333,  respectively.  Table  3  simply  shows  cross- 
sections  rather  than  efficiency  factors  for  turbule  radius  0.5  m  <  r  <  45  m,  wavelength 
X  =  15  m,  and  refractive  index  m  =  1.0000015.  The  ADF  formula  was  derived 
assuming  a  refractive  index  near  1  and  a  large  size  parameter  x  (that  is,  r  ►  X).  Table 
1  and  figure  3  show  clearly  good  agreement  between  program  ASCT  and  formula  results 
for  refractive  index  m  =  1.0000015.  For  refractive  index  m  =  1.33333,  program 
ASCT  and  formula  results  differ  slightly.  This  should  be  expected  since  the  refractive 
index  is  greater  than  1. 

Figure  4  shows  an  improvement  in  the  error  for  increasing  size  parameter  x.  This  is  in 
accordance  with  the  large  size  parameter  assumption  of  the  formula.  Using  the 
parameters  of  table  1  for  radius  of  0.5  m,  equation  (20)  gives  an  extinction  efficiency 
of  6.234E-13,  which  is  only  14.4  percent  different  than  the  program  results.  Using  the 
parameters  of  table  1  for  a  radius  of  0.05  m,  equation  (20)  gives  an  efficiency  of 
3.23320E-09  and  the  ASCT  program  gives  a  result  of  3.23328E-09,  indicating  the 
improved  agreement  for  smaller  size  parameter.  In  this  latter  case,  the  angular 
dependence  of  the  scattering  is  very  nearly  isotropic,  which  is  consistent  with  the  fact 
that  monopole  is  the  dominant  mode  of  scatter  when  x  is  small. 

The  phase  function  gives  out  information  about  the  angular  distribution  of  the  scattered 
wave.  Figure  5  shows  the  phase  function  curve  for  a  turbule  whose  radius  (r  =  30  m) 
is  large  compared  to  the  wavelength  (X  =  5  m)  for  refractive  indices  m  =  1.0000015 
and  m  =  1.33333.  Figure  6  shows  the  phase  function  curve  for  a  turbule  whose  radius 
(r  =  5  m)  is  small  compared  to  the  wavelength  (X  =  30  m).  Figure  7  shows  the  phase 
function  for  a  user-supplied  distribution  of  turbule  radii.  The  phase  function  for  a  log¬ 
normal  distribution  of  turbule  radii  is  shown  in  figure  8.  Tables  4  through  1 1  show  the 
corresponding  data  for  figures  5  through  8.  The  results  obtained  from  program  ASCT 
are  in  accord  with  those  predicted  by  the  "anomalous  diffraction"  formula  and  the  small 
size  approximation.  ASCT  correctly  models  the  acoustic  scattering  properties  of 
spherical  turbules. 

4.  Conclusion 

An  acoustic  scattering  code,  Program  ASCT,  is  now  available  for  determining  the 
scattering  properties  of  nonrelational  spherical  turbules.  It  can  generate  the  cross- 
sections  and  phase  function  for  four  modes  of  operation.  Program  ASCT  was  developed 
by  inserting  acoustic  interface  relations  in  the  electromagnetic  scattering  QGAUS  code. 
QGAUS  was,  in  turn,  developed  from  the  AGAUS  code  through  conversion  to  quadruple 
precision  arithmetic.  The  majority  of  the  modifications  to  program  QGAUS  were  in 
subroutine  MIEGX,  where  the  main  Mie  calculations  are  performed.  Modifications  to 
this  routine  were  based  on  the  derivations  provided  by  Dr.  G.  H.  Goedeeke.  These 
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derivations  were  compared  to  cross-section  and  phase  function  equations  given  in  Morse 
(1968)  for  scattering  by  spherical  particles.  Aside  from  different  notation  used,  the 
formulas  for  the  cross-sections  and  phase  function  were  found  to  be  identical.  The 
extinction  results  from  program  ASCT  were  then  compared  to  those  predicted  by  the 
"anomalous  diffraction"  formula  and  found  to  be  in  remarkably  good  accord.  Thus, 
there  is  now  available  an  accurate  code  for  determining  the  acoustic  scattering  properties 
of  spherical  turbules  modeled  as  uniform  acoustic  index-of-refraction  inhomogeneities. 
The  angular  dependence  for  small  size  parameter  was  compared  to  that  found  in  Morse 
(1968,  page  427).  Morse  shows  an  angular  dependence  with  predominant  scattering  in 
the  backward  direction.  ASCT  results  show  an  isotropic  scattering  pattern.  The 
difference  is  caused  by  the  different  boundary  conditions  for  the  two  derivations.  Morse 
considers  a  scatterer  with  a  definite  surface.  For  this  case,  the  boundary  conditions  are 
that  pressure  is  continuous  and  that  the  normal  velocity  is  continuous.  ASCT  considers 
a  scatterer  without  a  definite  surface.  In  the  ASCT  case,  the  boundary  conditions  are 
that  the  pressure  is  continuous  and  that  the  mass  flow  (density  times  velocity)  is 
continuous.  Morse  shows  the  dipole  term  contributes  to  the  lowest  order  of  scattering. 
In  ASCT,  only  the  monopole  term  contributes  to  the  lowest  order  of  scattering. 
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Figure  3.  Percent  error  in  ASCT  and  ADF  formula  results  for  refractive 
index  m  =  1.0000015. 
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Figure  5.  Phase  function  for  a  turbule  whose  radius  is  small  compared 
to  wavelength. 


Figure  6.  Phase  function  for  a  turbule  whose  radius  is  small  compared 
to  wavelength. 
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Figure  7.  Phase  function  for  arbitray  user-supplied  model. 


Figure  8.  Phase  function  for  tog-normal  size  distribution  model. 
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TABLE  1.  TABLE  GENERATION  MODE 


m  -  1.0000015 


RLO  -  5 . 00000E-01 .  RHI  -  1.05000E+02.  DELR  -  5.00000E+00 

Wavelength  -  5.00  m 

Real  Index  -  1.00000150  Imag.  Indjx  -  .00000000 


Size 

_ SXslI 

Parameter 

1— .  HU.fi 

W  <11 > 1 1 

)  ADF 

.50000  6 . 28319E-01  5.33358E-13  5.33358E-13  0.00000E+00  1.77653E-12 

5.50000  6 . 91150E+00  2.13878E-10  2.13878E-10  0.00000E+00  2 . 14960E-10 

10.50000  1 . 31947E+01  7.82349E-10  7.82349E-10  0.00000E+00  7.83449E-10 

15.50000  1.94779E+01  1.70614E-09  1.70614E-09  O.OOOOOE+OO  1.70724E-09 

20.50000  2 . 57611E+01  2.98524E-09  2.98524E-09  O.OOOOOE+OO  2.98634E-09 

25.50000  3 . 20442E+01  4.61964E-09  4.61964E-09  O.OOOOOE+OO  4.62075E-09 

30.50000  3 . 83274E+01  6.60936E-09  6.60936E-09  O.OOOOOE+OO  6.61046E-09 

35.50000  4.46106E+01  8.95435E-09  8.95435E-09  O.OOOOOE+OO  8.95548E-09 

40.50000  5.08938E+01  1.16547E-08  1.16547E-08  O.OOOOOE+OO  1.16558E-08 

45.50000  5 . 71770E+01  1.47103E-08  1.47103E-08  O.OOOOOE+OO  1.47114E-08 

50.50000  6 . 34602E+01  1.81212E-08  1.81212E-08  O.OOOOOE+OO  1.81224E-08 

55.50000  6 . 97434E+01  2.18873E-08  2.18873E-08  O.OOOOOE+OO  2.18886E-08 

60.50000  7 . 60265E+01  2.60089E-08  2.60089E-08  O.OOOOOE+OO  2.60102E-08 

65.50000  8 . 23097E+01  3.04858E-08  3.04858E-08  O.OOOOOE+OO  3.04870E-08 

70.50000  8 . 85929E+01  3.53176E-08  3.53176E-08  O.OOOOOE+OO  3.53192E-08 

75.50000  9.48761E+01  4.05039E-08  4.05039E-08  O.OOOOOE+OO  4.05066E-08 

80.50000  1 . 01159E+02  4.60471E-08  4.60471E-08  O.OOOOOE+OO  4.60494E-08 

85.50000  1 . 07442E+02  5.19455E-08  5.19455E-08  O.OOOOOE+OO  5.19475E-08 

90.50000  1 . 13726E+02  5.81974E-08  5.81974E-08  O.OOOOOE+OO  5.82009E-08 

95.50000  1 . 20009E+02  6.48068E-08  6.48068E-08  O.OOOOOE+OO  6.48095E-08 

100.50000  1.26292E+02  7.17713E-08  7.17713E-08  O.OOOOOE+OO  7.17735E-08 
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TABLE  2.  TABLE  GENERATION  MODE 


m  -  1.33333 


RLO  -  5 . 00000E-01 .  RHI  -  1.05000E+02.  DELR  -  5.00000E+00 

Wavelength  -  5.00  m 

Real  Index  -  1.33333000  Imag-  Index  -  .00000000 


Size 

_ RCnO  Parameter  Q(EXT)  Q(SCA)  Q(ABS)  .  APF  . 

.50000  6 . 28319E-01  4.44193E-02  4.44193E-02  0.00000E+00  8.68767E-02 

5.50000  6 . 91150E+00  3.81237E+00  3.81237E+00  O.OOOOOE+OO  3.07148E+00 

10.50000  1. 31947E+01  2.33832E+00  2.33832E+00  O.OOOOOE+OO  1.82620E+00 

15.50000  1 . 94779E+01  2.00868E+00  2.00868E+00  0.00000E+00  1.87679E+00 

20.50000  2 . 57611E+01  2.56390E+00  2.56390E+00  0.00000E+00  2.24661E+00 

25.50000  3 . 20442E+01  2.07078E+00  2.07078E+00  O.OOOOOE+OO  1.90576E+00 

30.50000  3 . 83274E+01  2.05790E+00  2.05790E+00  O.OOOOOE+OO  1.93689E+00 

35.50000  4.46106E+01  2.33004E+00  2.33004E+00  O.OOOOOE+OO  2.13875E+00 

40.50000  5 . 08938E+01  2.14009E+00  2.14009E+00  O.OOOOOE+OO  1.93696E+00 

45.50000  5 . 71770E+01  2.08112E+00  2.08112E+00  O.OOOOOE+OO  1.95759E+00 

50.50000  6 . 34602E+01  2.24026E+00  2.24026E+00  O.OOOOOE+OO  2.09650E+00 

55.50000  6 . 97434E+01  2.05988E+00  2.05988E+00  O.OOOOOE+OO  1.95275E+00 

60.50000  7 . 60265E+01  2.06883E+00  2.06883E+00  O.OOOOOE+OO  1.96807E+00 

65.50000  8 . 23097E+01  2.19703E+00  2.19703E+00  O.OOOOOE+OO  2.07396E+00 

70.50000  8 . 85929E+01  2.04687E+00  2.04687E+00  O.OOOOOE+OO  1.96223E+00 

75.50000  9.48761E+01  2 . 11949E+00  2.11949E+00  O.OOOOOE+OO  1.97440E+00 

80.50000  1 . 01159E+02  2.16265E+00  2.16265E+00  O.OOOOOE+OO  2.05996E+00 

85.50000  1.07442E+02  2.04734E+00  2.04734E+00  O.OOOOOE+OO  1.96855E+00 

90.50000  1 . 13726E+02  2.05237E+00  2.05237E+00  O.OOOOOE+OO  1.97864E+00 

95.50000  1.20009E+02  2.14777E+00  2.14777E+00  O.OOOOOE+OO  2.05041E+00 

100.50000  1 . 26292E+02  2.04060E+00  2.04060E+00  O.OOOOOE+OO  1.97306E+00 


22 


TABLE  3.  TABLE  GENERATION  MODE  (CROSS-SECTIONS) 
m  -  1.0000015 


RLO  -  5.00000E-01.  RHI  -  4.50000E+01.  DELR  -  1. OOOOOE+OO 

Wavelength  —  15.00  m 

Real  Index  -  1.00000150  Imag.  Index  -  .00000000 


R(m) 


Size 

Parameter 


Cross  Sections  In  Square  Meters 


Q(EXT) 


Q(SCA) 


Q(AB.S), 


ADF 


.50000  2 . 09440E-01 

1.50000  6 . 28319E-01 

2.50000  1 . 04720E+00 

3.50000  1 . 46608E+00 

4.50000  1 . 88496E+00 

5.50000  2 . 30383E+00 

6.50000  2 . 72271E+00 

7.50000  3 . 14159E+00 

8 . 50000  3 . 56047E+00 

9.50000  3 . 97935E+00 

10.50000  4 . 39823E+00 

11.50000  4 . 81711E+00 

12.50000  5 . 23599E+00 

13.50000  5 . 65487E+00 

14.50000  6.07375E+00 

15.50000  6.49262E+00 

16.50000  6 . 91150E+00 

17.50000  7 . 33038E+00 

18.50000  7 . 74926E+00 

19.50000  8 . 16814E+00 

20.50000  8 . 58702E+00 

21.50000  9 . 00590E+00 

22.50000  9 . 42478E+00 

23.50000  9 . 84366E+00 

24.50000  1 . 02625E+01 

25.50000  1 . 06814E+01 

26.50000  1 . 11003E+01 

27.50000  1 . 15192E+01 

28.50000  1 . 19381E+01 

29.50000  1 . 23569E+01 

30.50000  1 . 27758E+01 

31.50000  1 . 31947E+01 

32.50000  1 . 36136E+01 

33.50000  1 . 40324E+01 

34.50000  1 .44513E+01 

35.50000  1.48702E+01 

36.50000  1 . 52891E+01 

37.50000  1 . 57080E+01 

38.50000  1 . 61268E+01 

39.50000  1 . 65457E+01 

40.50000  1 . 69646E+01 


5 . 93984E-15 
3.77008E-12 
6 . 18948E-11 
3 . 22713E-10 
9.61909E-10 
2.16275E-09 
4. 24848E-09 
7 . 64968E-09 
1 . 27259E-08 
1 . 98716E-08 
2 . 97172E-08 
4 . 29363E-08 
6 . 00489E-08 
8 . 17080E-08 
1.08862E-07 
1 . 42374E-07 
1 . 82929E-07 
2.31494E-07 
2 . 89313E-07 
3 . 57389E-07 
4. 36610E-07 
5 . 28296E-07 
6 . 33937E-07 
7 . 54649E-07 
8 . 91579E-07 
1 . 04642E-06 
1 . 22084E-06 
1 . 41606E-06 
1.63357E-06 
1 . 87541E-06 
2 . 14335E-06 
2.43877E-06 
2 . 76356E-06 
3 . 12004E-06 
3 . 51005E-06 
3 . 93519E-06 
4 . 39776E-06 
4 . 90033E-06 
5.44475E-06 
6.03294E-06 
6 . 66764E-06 


5 . 93984E-15 
3 . 77008E-12 
6 . 18948E-11 
3 . 22713E-10 
9 . 61909E-1Q 
2 . 16275E-09 
4 . 24848E-09 
7 . 64968E-09 
1. 27259E-08 
1 . 98716E-08 
2 . 97172E-08 
4. 29363E-08 
6 . 00489E-08 
8 . 17080E-08 
1.08862E-07 
1.42374E-07 
1.82929E-07 
2.31494E-07 
2 . 89313E-07 
3 . 57389E-07 
4. 36610E-07 
5 . 28296E-07 
6 . 33937E-07 
7 . 54649E-07 
8 . 91579E-07 
1 . 04642E-06 
1 . 22084E-06 
1 , 41606E-06 
1 . 63357E-06 
1 . 87541E-06 
2 . 14335E-06 

2.76356E-06 
3 . 12004E-06 
3 . 51005E-06 
3 . 93519E-06 
4 . 39776E-06 
4 . 90033E-06 
5.44475E-06 
6 . 03294E-06 
6 . 66764E-06 


0 . Q0000E+00 
0 . 00000E+00 
0 . OOOOOE+OO 
0 . 00000E+00 
0 . 00000E+00 
0. 00000E+00 
0 . 00000E+00 
0 . 00000E+00 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0.  OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 . OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 


1 . 55031E-13 
1 . 25575E-11 
9 . 68946E-11 
3 . 72230E-10 
1 .01716E-09 
2 . 26981E-09 
4.42785E-09 
7 . 84846E-09 
1 . 29484E-08 
2 . 02038E-08 
3.01507E-08 
4. 33841E— 08 
6 . 05591E-08 
8 . 23900E-08 
1 . 09651E— 07 
1 . 43175E— 07 
1 . 83855E-07 
2 . 32644E— 07 
2 . 90554E— 07 
3 . 58656E-07 
4 . 38082E-07 
5 . 30021E-07 
6 . 35726E-07 
7 . 56504E-07 
8 . 93725E-07 
1 . 04882E-06 
1 . 22327E-06 
1.41863E-06 
1 . 63651E-06 
1 . 87857E-06 
2 . 14654E-06 
2.44220E-06 
2 . 76741E— 06 
3 . 12406E-06 
3 . 51412E-06 
3 . 93961E-06 
4 . 40262E-06 
4 . 90529E-06 
5 . 44982E-06 
6 . 03848E-06 
6 . 67359E-06 
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TABLE  4.  SINGLE  TURBULE  MODEL 


A  <  7 

m  -  1.0000015 


Wavelength 
Real  Index 

-  5.00  m 

-  1.00000150  Imag.-  Index  -  .00000000 

R(m) 

Size 

Parameter  O(EXT)  O(SCA)  O(ABS) 

ADF 

30.00000  3 

. 76991E+01  6 . 39439E-09  6.39439E-09  0.00000E+00 

6 . 39550E-09 

***  Phase  Function  (Integral  over  solid  angle  is  unity)  *** 


Angle 

(Dee). 

Phase 

Function 

Angle 

(Dez) 

Phase 

Function 

Angle 

(Dez) 

Phase 

Function 

.000 

1.00549E+02 

.500 

9 . 83923E+01 

1.000 

9 . 21591E+01 

2.000 

7 . 04779E+01 

3.000 

4. 39194E+01 

4.000 

2 . 11377E+01 

5.000 

6 . 90125E+00 

6.000 

9 . 74323E-01 

8.000 

5 . 51044E— 01 

10.000 

4.06732E-01 

12.000 

5 . 57109E-03 

14.000 

1 . 26362E— 01 

16.000 

1 . 18102E-02 

18.000 

2 . 81551E-02 

20.000 

2 . 10238E-02 

24 . 000 

1 . 50308E-02 

28.000 

5 . 92785E-03 

32.000 

7 . 73665E-04 

36.000 

1.46115E-04 

40.000 

1 . 21146E-03 

44 . 000 

1.42367E-03 

48.000 

5 . 82885E-04 

52.000 

6 . 96884E-06 

56.000 

2.40956E-04 

60.000 

4.48025E-04 

64.000 

1 . 51021E-04 

68.000 

1 . 43342E-05 

72.000 

2.05845E-04 

76.000 

1 . 17420E-04 

80.000 

7 . 05946E-06 

84.000 

1 . 33905E-04 

88.000 

3 . 38229E-05 

92.000 

4 . 57344E-05 

96.000 

7.08101E-05 

100.000 

9 . 23468E-06 

104.000 

6 . 78840E-05 

108.000 

3 . 88159E-06 

112.000 

5 . 38492E-05 

116.000 

1 . 00587E-05 

120.000 

3 . 10136E-05 

124.000 

3 . 07702E-05 

128.000 

2 . 38710E-06 

132.000 

3 . 82840E-05 

136.000 

1 . 81152E-05 

140.000 

1 . 14544E-06 

142.000 

1 . 17876E-05 

144.000 

2 . 53459E-05 

146.000 

3 . 28441E— 05 

148.000 

3 . 10271E-05 

150.000 

2 . 23238E-05 

152.000 

1 . 17272E-05 

154.000 

3 . 61213E-06 

156.000 

1 . 21887E-07 

158.000 

1 . 18045E-06 

160.000 

5 . 37666E-06 

162.000 

1 . 09503E-05 

164.000 

1 . 64533E-05 

166.000 

2 . 10105E-05 

168.000 

2.42952E-05 

170.000 

2 . 63621E-05 

172.000 

2.74691E-05 

174.000 

2 . 79327E-05 

176.000 

2 . 80440E-05 

178.000 

2 . 80226E-05 

180.000 

2 . 80014E-05 

Test  integral  of  phase  function  is 

:  1 . 04217E+00 

TABLE  5.  SINGLE  TORBULE  MODEL 


A  <  7 

m  -  1,33333 


Wavelength  -  5.00  m 

Real  Index  -  1,33333000  Imag.  Index  -  .00000000 


Size 

Rfm)  Parameter 

0(EXT) 

Of'SCA'l 

Of ABS^ 

ADF 

30.00000  3 . 76991E+01 

2 . 13171E+00 

2 . 13171E+00 

0 . 00000E+00 

2.00004E+00 

***  Phase  Function  (Integral  over  solid  angle  is  unity)  *** 


Angle 

(Deg) 

.000 

6 . 19061E+01 

2.000 

3 . 79818E+01 

5.000 

6 . 25238E— 01 

10.000 

8 . 24598E— 01 

16.000 

6 . 80251E— 01 

24.000 

1 . 99759E-01 

36.000 

3 . 36434E— 02 

48.000 

7.44567E-02 

60.000 

2 . 60021E-02 

72.000 

2 . 76854E-03 

84.000 

1 . 00182E-02 

96.000 

4. 77071E-03 

108.000 

9 . 80708E-03 

120.000 

1 . 79242E-04 

132.000 

3 . 67086E-03 

142.000 

3 . 54139E-02 

148.000 

3 . 13271E-02 

154.000 

8 . 61732E-03 

160.000 

1 . 25980E-02 

166.000 

3 . 12651E-02 

172.000 

4 . 45448E-03 

178.000 

3 . 72092E-02 

Angle 

Phase 

iPH 

.500  6 . 01188E+01 


3.000 

1 . 91703E+01 

6.000 

6 . 73745E-01 

12.000 

3 . 15137E-02 

18.000 

3 . 13066E-01 

28.000 

3 . 18359E-01 

40.000 

5 . 58591E-02 

52.000 

2 . 13933E-02 

64.000 

3 . 11628E-02 

76.000 

5 . 94687E-03 

88 . 000 

1 . 30953E-03 

100 . 000 

1 . 03444E-03 

112.000 

9 . 99972E-03 

124.000 

1 . 85665E-03 

136.000 

7 . 12088E-03 

144.000 

4 . 36453E-02 

150.000 

2 . 95026E-02 

156.000 

1 . 33666E-03 

162.000 

1 . 12088E-02 

168.000 

3 . 94053E-03 

174.000 

6 . 24914E-02 

180.000 

1 .46867E-01 

Angle 

(Dez) 

Phase 

Function 

1.000 

5 . 50098E+01 

4.000 

5 . 94872E+00 

8.000 

2 . 49006E+00 

14 . 000 

2 . 42282E-01 

20.000 

2 . 85330E-01 

32.000 

1 . 65788E-01 

44.000 

1 . 04645E-01 

56.000 

4 . 99427E-03 

68.000 

1 . 57804E— 02 

80.000 

1.43038E-02 

92.000 

2 . 38238E-03 

104.000 

1 . 82243E-03 

116.000 

2 . 68559E-03 

128.000 

1 . 07126E-02 

140.000 

1 . 08049E-02 

146.000 

2 . 30289E-02 

152.000 

3 . 07851E-02 

158.000 

3 . 37359E-03 

164.000 

1 . 77247E-02 

170.000 

4 . 49724E-03 

176.000 

2 . 87252E-02 

Test  integral  of  phase  function  is:  1.00015E+00 
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TABLE  6.  SINGLE  TURBULE  MODEL 


A  <  7 

m  -  1.0000015 


Wavelength  -  30.00  m 

Real  Index  -  1.00000150  Imag.  Index  -  .0C000000 


R(m) 

Size 

Parameter 

Q(EXT) 

0(SCA) 

0(ABS> 

ADF 

5 . 00000 

1.04720E+00 

3.15228E-12 

3 . 15228E-12 

0 . 00000E+00 

4 . 93480F.--12 

***  Phase  Function  (Integral  over  solid  angle  is  unity)  *** 


Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

(Dee) 

Function 

(Dee) 

Function 

(Dee) 

Function 

.000 

1 . 21434E— 01 

.500 

1 . 21432E-01 

1.000 

1 . 21426E— 01 

2.000 

1 . 21402E-01 

3.000 

1 . 21361E-01 

4.000 

1 . 21305E— 01 

5.000 

1 . 21232E-01 

6.000 

1 . 21143E-01 

8.000 

1.20917E-01 

10.000 

1 . 20627E-01 

12.000 

1 . 20275E-01 

14.000 

1 . 19861E-01 

16.000 

1 . 19386E-01 

18.000 

1 . 18851E-01 

20.000 

1.18258E— 01 

24.000 

1 . 16903E— 01 

28.000 

1 . 15335E-01 

32.000 

1 . 13568E— 01 

36.000 

1.11619E— 01 

40.000 

1.09507E-01 

44.000 

1 . 07250E-01 

48.000 

1 . 04868E-01 

52.000 

1 . 02381E-01 

56.000 

9 . 98096E-02 

60.000 

9 . 71734E-02 

64.000 

9.44920E-02 

68.000 

9 . 17842E— 02 

72.000 

8 . 90681E-02 

76.000 

3 . 63606E-02 

80.000 

8 . 36776E-02 

84.000 

8 . 10337E-02 

88.000 

7 . 84423E-02 

92.000 

7 . 59153E-02 

96.000 

7 . 34634E-02 

100 . 000 

7 . 10959E-02 

104.000 

6 . 88210E— 02 

108.000 

6 . 66455E-02 

112.000 

6 . 45751E-02 

115.000 

6 . 26145E-02 

120.000 

6 . 07673E-02 

124.000 

5 . 90365E-02 

128.000 

5 . 74241E-02 

132.000 

5 . 59315E-02 

136.000 

5 . 45596E-02 

140.000 

5 . 33089E-02 

142.000 

5 . 27291E-02 

144.000 

5 . 21795E-02 

146 . 000 

5 . 16601E-02 

148.000 

5 . 11710E-02 

150.000 

5 . 07120E-02 

152.000 

5 . 02831E-02 

154.000 

4 . 98842E— 02 

156.000 

4 . 95153E-02 

158.000 

4 . 91762E-02 

160.000 

4 . 88669E-02 

162.000 

4 . 85873E-02 

164.000 

4 . 83374E-02 

166.000 

4 . 81170E-02 

168.000 

4. 79262E-02 

170.000 

4 . 77648E-02 

172.000 

4 . 76329E-02 

174.000 

4 . 75303E-02 

176.000 

4 . 74570E-02 

178.000 

4.74131E— 02 

180.000 

4 . 73984E-02 

Test  integral  of  phase  function  is:  1.00005E+00 
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TABLE  7.  SINGLE  TURBULE  MODEL 


X  <  7 

m  -  1.33333 


Wavelength  -  30.00  m 

Real  Index  -  1.33333000  Imag.  Index  -  .00000000 


Size 

R(m)  Parameter 

ADF 

5.00000  1 . 04720E+00 

3 . 08358E-01 

3 . 08358E-01  O.OOOOOE+OO  2 

.  37171E-01 

***  Phase  Function  (Integral 

over  solid  angle  is  unity) 

Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

WISH 

Function 

SnNuTaHip 

Function 

.000 

1 . 16616E-01 

.500 

1 . 16615E-01 

1.000 

1 . 16609E— 01 

2.000 

1 . 16588E— 01 

3.000 

1 . 16553E-01 

4.000 

1 . 16503E-01 

5.000 

1 . 16439E— 01 

6.000 

1 . 16361E-01 

8.000 

1. 16164E-01 

10.000 

1 . 15910E-01 

12  noo 

1 . 15602E-01 

14 . 000 

1 . 15240E-01 

16.000 

1 . 14824E— 01 

18.0.0 

1 . 14356E-01 

20.000 

1 . 13836E-01 

24.000 

1 . 12650E-01 

28.000 

1.11275E-01 

32.000 

1 . 09725E-01 

36 . 000 

1 . 08015E-01 

40 . 000 

1 . 06160E-01 

44 . 000 

1 . 04176E-01 

48.000 

1 . 02081E-01 

52.000 

9 . 98906E-02 

56.000 

9 . 76235E-02 

60.000 

9 . 52967E-02 

64.000 

9 . 29272E-02 

68.000 

9 . 05313E-02 

72.000 

8 . 81250E-02 

76.000 

8 . 57231E-02 

80.000 

8 , 33397E-02 

84.000 

8 . 09876E-02 

88.000 

7 . 86789E-02 

92.000 

7 . 64244E-02 

96.000 

7 . 42336E-02 

100.000 

7 . 21151E-02 

104.000 

7 . 00765E-02 

108 . 000 

6 . 81240E-02 

112.000 

6 , 62631E-02 

116.000 

6 . 44984E-02 

120.000 

6 . 28336E-02 

124.000 

6 . 12714E-02 

128.000 

5 . 98142E-02 

132.000 

5 . 84636E-02 

136.000 

5 . 72208E-02 

140.000 

5 . 60865E-02 

142.000 

5 . 55601E-02 

144.000 

5 . 50610E-02 

146.000 

5 . 45891E-02 

148.000 

5 . 41445E-02 

150.000 

5 . 37271E-02 

152.000 

5 . 33369E-02 

154.000 

5 . 29739E-02 

156.000 

5 . 26379E-02 

158.000 

5 . 23291E-02 

160.000 

5.20473E-02 

162.000 

5 . 17926E-02 

164.000 

5 . 15647E-02 

166.000 

5 . 13638E-02 

168.000 

5 . 11898E-02 

170.000 

5 . 10426E-02 

172.000 

5 . 09222E-02 

174.000 

5 . 08286E-02 

176.000 

5 . 07618E-02 

178,000 

5.07217E-02 

180.000 

5 . 07083E-02 

Test  integral  of  phase  function  is:  1.00Q04E+00 
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TABLE  8.  ARBITRARY  OSER-SUPPLIED  MODEL 
m  -  1.0000015 


Radius  Relative  No,  Radius  Relative  No.  Rajivs _  Relative  No. 

1.900000E-01  1 . 000000E+00  2.500000E-01  1.000000E+00  3.000000E-01  1.000000E+00 
3 . 500000E-01  1.000000E+00  4.000000E-01  1.000000E+00 

Average  numerical  dry  volume  per  particle  is:  1.211598E-01  cubic  meters 


Wavelength  -  .11m 

Real  Index  -  1.00000150 

Imag .  Index  -  . 00000000 

Number  Density  -  2.10000E-01  part/cubic  meter 


Size 

Maximum 

GrouD  cross 

sections  in  sauare  meters 

Der  oarticle 

Group 

Radius  Extinction 

Absorption 

Scattering  Backscatter 

1 

.400  4. 38812E-10 

4.48416E-44 

4 . 38812E— 

10  3 . 18346E-13 

***  Phase  Function  (integral 

over  solid  angle  is  unity)  *** 

Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

(Dee) 

Function 

(Dee) 

Function 

(Dee) 

Function 

.000 

2.69153E+01 

2.813 

2 . 20677E+01 

5.625 

1 . 18638E+01 

8.438 

3 . 84435E+00 

11.250 

6 . 91398E-01 

14.063 

2 . 09543E-01 

16.875 

1 . 54649E-01 

19.688 

6 . 20410E-02 

22.500 

3 . 56652E-02 

25.313 

3.06032E-02 

28.125 

1 . 63331E-02 

30.938 

1 . 10230E-02 

33.750 

9 . 18239E-03 

36.563 

5.49717E-03 

39.375 

4 . 27419E-03 

42.188 

3 . 85865E-03 

45.000 

2 . 64059E-03 

47.813 

2 . 11049E-03 

50.625 

1 . 89161E-03 

53.438 

1 . 36339E-03 

56.250 

1 . 10207E-03 

59.063 

1 . 04922E-03 

61.875 

8. 42690E-04 

64.688 

6 . 78473E-04 

67.500 

6 . 40439E-04 

70.313 

5 . 48956E-04 

73.125 

4 . 33446E-04 

75.938 

3 . 97018E-04 

78.750 

3 . 76826E-04 

81.563 

3 . 19339E-04 

84.375 

2 . 75281E-04 

87.188 

2 . 63894E-04 

90.000 

2.44566E-04 

92.813 

2 . 08157E-04 

95.625 

1 . 83319E-04 

98.438 

1 . 76756E-04 

101.250 

1 . 68788E-04 

104.063 

1 . 51246E-04 

106.875 

1 . 35169E-04 

109.688 

1 . 28484E-04 

112.500 

1 . 25767E-04 

115.313 

1 . 19032E-04 

118.125 

1 . 08046E-04 

120.938 

9 . 81422E-05 

123.750 

9 . 28082E-05 

126.563 

9 . 08515E-05 

129.375 

8 . 90175E-05 

132.188 

8 . 53915E-05 

135.000 

8 . 03582E-05 

137.813 

7.54421E-05 

140.625 

7 . 18429E-05 

143.438 

6 . 98105E-05 

146.250 

6 . 88598E-05 

149.063 

6 . 82929E-05 

151.875 

6 . 75910E-05 

154.688 

6 . 65425E-05 

157.500 

6 . 51827E-05 

160.313 

6 . 36674E-05 

163.125 

6 . 21686E-05 

165.938 

6 . 08176E-05 

168.750 

5 . 96912E-05 

171.563 

5 . 88209E-05 

174.375 

5 . 82100E-05 

177.188 

5 . 78500E-05 

180.000 

5 . 77313E-05 

Test  integral  of  phase  function  is:  1.06681E+00 
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TABLE  9.  ARBITRARY  USER- SUPPLIED  MODEL 


m  -  1.33333 

_ Radju? _  fasUtjvg  HSU  Radius  Relative  No.  Radius  Relative  No, 

1 . 900000E-01  1 . 000000E+00  2.500000E-01  1.000000E+00  3.000000E-01  1.000000E+00 
3 . 500000E-01  1 . 000000E+00  4.000000E-01  1.000000E+00 

Average  numerical  dry  volume  per  particle  is:  1.211598E-01  cubic  meters 


Wavelength 
Real  Index 
Imag .  Index 
Number  Density 


.11  m 

1.33333000 
. 00000000 

2.10000E-01  part/cubic  meter 


Size 

Maximum 

Grouo  cross 

sections  in  sauare  meters 

Der  particle 

Group 

Radius  Extinction 

Absorption 

Scattering  Backscatter 

1 

.400  6.42572E-01 

0 . OOOOOE+OO 

6 . 42572E- 

01  1 . 76685E+00 

***  Phase  Function  (integral 

over  solid  angle  is  unity)  *** 

Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

_(Dee) 

Function 

(Dee) 

Function 

(Dee) 

Function 

.000 

1 . 56191E+01 

2.813 

1 . 19070E+01 

5.625 

5 . 20701E+00 

8.438 

1 . 55340E+00 

11.250 

8 . 72676E-01 

14.063 

6 . 94558E-01 

16.875 

4 . 63789E-01 

19.688 

3 . 92105E-01 

22.500 

3 . 25768E-01 

25.313 

2 . 19202E-01 

28.125 

1 . 78952E-01 

30.938 

1 . 68788E-01 

33.750 

1 . 32202E-01 

36.563 

1 . 08059E-01 

39.375 

9 . 91544E-02 

42.188 

7 . 95723E-02 

45.000 

6.49130E-02 

47.813 

5 . 96756E-02 

50.625 

4 . 94084E-02 

53.438 

4 . 11315E-02 

56.250 

3 . 74403E-02 

59.063 

2 . 99048E-02 

61.875 

2.49234E-02 

64.688 

2 . 31956E-02 

67.500 

1 . 85913E-02 

70.313 

1 . 60056E-02 

73.125 

1 . 50842E— 02 

75.938 

1 . 17547E-02 

78.750 

9 . 95956E-03 

81.563 

9 . 55965E-03 

84.375 

7 . 62468E-03 

87.188 

6 . 84997E-03 

90 . 000 

6 . 63866E-03 

92.813 

5.49053E-03 

95.625 

5 . 30885E-03 

98.438 

5 . 31182E-03 

101.250 

5 . 21668E-03 

104.063 

5.40052E-03 

106.875 

5 . 47881E-03 

109.688 

5 . 56325E-03 

112.500 

5.45785E-03 

115.313 

5 . 58924E-03 

118.125 

6 . 15020E-03 

120.938 

6.34299E-03 

123.750 

6 . 352S8E-03 

126.563 

7 . 99460E-03 

129.375 

9 . 02378E-03 

132.188 

1 . 10437E-02 

135.000 

1 . 53470E-02 

137.813 

1 . 62270E-02 

140.625 

2 . 27691E-02 

143.438 

2 . 94942E-02 

146.250 

2 . 78958E-02 

149.063 

3 . 62969E-02 

151.875 

3 . 89313E-02 

154.688 

2 . 87069E-02 

157.500 

2 . 53064E-02 

160.313 

1 . 94413E-02 

163.125 

1 . 45743E-02 

165.938 

1 . 49503E-02 

168.750 

2 . 14306E-02 

171.563 

2 . 16007E-02 

174.375 

3 . 05357E-02 

177.188 

1 . 32714E-01 

180.000 

2 . 18811E-01 

Test  integral  of  phase  function  is:  1.03898E+00 
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TABLE  10.  LOG-NORMAL  SIZE  DISTRIBUTION  MODEL 
m  -  1.0000015 

RBAR  -  2.25000E-01.  SIGMA  -  8.00000E-02.  RLO  -  2.00000E-01 
RHI  -  6.50000E-01.  DENS  -  5.00000E+00 

Average  numerical  dry  volume  per  particle  is:  3.086332E-01  cubic  meters 


Wavelength  - 
Real  Index  - 
Imag .  Index  - 
Number  Density  - 


.12  m 

1.00000150 

.00000000 

5.00000E+00  part/cubic  meter 


Size 

Maximum 

Grouo  cross 

sections  in  i 

sauare  meters 

per  particle 

GrouD 

Radius  Extinction 

Absorption 

Scattering  Backscatter 

1 

.650  1.41035E-09 

0 .  OOOOOE+OO 

1.41035E- 

09  5 . 63142E-13 

***  Phase  Function  (Integral 

over  solid  angle  is  unity)  *** 

Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

(Deg) 

Function 

(Deg) 

Function 

(Dee) 

Function 

.000 

5.46813E+01 

2.813 

3 . 57618E+01 

5.625 

9 . 18735E+00 

8.438 

1.1I536E+00 

11.250 

4. 41477E-01 

14.063 

1 . 31935E-01 

16.875 

7 . 58641E-02 

19.688 

3 . 32985E-02 

22.500 

2 . 35993E-02 

25.313 

1 . 31178E-02 

28.125 

9.94893E-03 

30.938 

6 . 10393E-03 

33.750 

4. 71830E-03 

36.563 

3 . 12050E-03 

39.375 

2 . 54196E-03 

42.188 

1 . 82794E— 03 

45.000 

1. 55004E-03 

47.813 

1.15871E— 03 

50.625 

9 . 90197E-04 

53.438 

7 . 55367E-04 

56.250 

6 . 57504E-04 

59.063 

5 . 22778E-04 

61.875 

4. 64370E-04 

64.688 

3 . 85438E-04 

67.500 

3 . 38978E-04 

70.313 

2.91634E-04 

73.125 

2 . 50455E-04 

75.938 

2 . 24592E-04 

78.750 

1 . 92143E-04 

81.563 

1 . 76692E-04 

84.375 

1 . 55715E-04 

87.188 

1. 40705 E-04 

90.000 

1 . 29404E-04 

92.813 

1 . 14679E-04 

95.625 

1 . 06163E-04 

98.438 

9 . 71798E-05 

101.250 

8 . 78992E-05 

104.063 

8 . 25674E-05 

106.875 

7 . 66793E-05 

109.688 

7.05860E-05 

112.500 

6 . 67580E-05 

115.313 

6 . 30485E-05 

118.125 

5 . 86993E-05 

120.938 

5 . 53438E-05 

123.750 

5 . 28847E-05 

126.563 

5 . 02146E-05 

129.375 

4. 74549E-05 

132.188 

4 . 52723E-05 

135.000 

4 . 36412E-05 

137.813 

4 . 21141E— 05 

140.625 

4 . 05212E-05 

143.438 

3 . 90216E-05 

146.250 

3 . 77783E-05 

149.063 

3 . 67912E-05 

151.875 

3 . 59630E-05 

154.688 

3 . 52097E-05 

157.500 

3 . 45032E-05 

160.313 

3 . 38544E-05 

163.125 

3 . 32842E-05 

165.938 

3 , 28063E-05 

168.750 

3 . 24241E-05 

171.563 

3 . 21348E-05 

174.375 

3 . 19330E-05 

177.188 

3 . 18140E-05 

180.000 

3 . 17748E-05 

Test  integral  of  phase  function  is:  1.13438E+00 


TABLE  11.  LOG-NORMAL  SIZE  DISTRIBUTION  MODEL 
m  -  1.33333 

RBAR  -  2 . 25000E-01.  SIGMA  -  8.00000E-02.  RLO  -  2.00000E-01 
RHI  -  6 . 50000E-01.  DENS  -  5.00000E+00 

Average  numerical  dry  volume  per  particle  is:  3.086332E-01  cubic  meters 


Wavelength  - 
Real  Index  - 
Imag .  Index  - 
Number  Density  - 


.12  m 

1.33333000 

.00000000 

5.00000E+00  part/cubic  meter 


Size  Maximum  Group  cross  sections  in  square  meters  per  particle 

Group  Radius  Extinction  Absorption  Scattering  Backscatter 

1  .650  1 . 11976E+00  O.OOOOOE+OO  1.11976E+00  2.59044E+00 


***  Phase  Function  (Integral  over  solid  angle  is  unity)  *** 


Angle 

Phase 

Angle 

Phase 

Angle 

Phase 

(Dee) 

Function 

(Dee) 

Function 

(Dee) 

Function 

.000 

2 . 79135E+01 

2.813 

1 . 68235E+01 

5.625 

3 . 92475E+00 

8.438 

1 . 28005E+00 

11.250 

7 . 36541E-01 

14.063 

5 . 85671E— 01 

16.875 

4 . 29428E-01 

19.688 

3.47377E-01 

22.500 

2 . 89255E— 01 

25.313 

2 . 35830E-01 

28.125 

2 . 02544E-01 

30.938 

1 . 61005E-01 

33.750 

1.42696E— 01 

36.563 

1 . 16183E-01 

39.375 

9 . 90076E-02 

42.188 

7 . 84624E-02 

45.000 

6 . 84321E-02 

47.813 

5 . 42196E-02 

50.625 

4. 72283E-02 

53.438 

3 . 81453E-02 

56.250 

3 . 31485E-02 

59.063 

2 . 71950E-02 

61.875 

2 . 27121E-02 

64.688 

1 . 86862E-02 

67.500 

1 . 62634E-02 

70.313 

1.33752E-02 

73.125 

1 . 15372E-02 

75.938 

9 . 96603E-03 

78.750 

8.47404E-03 

81.563 

7 . 61079E-03 

84.375 

6.65180E— 03 

87.188 

5 . 57197E-03 

90.000 

5.58129E-03 

92.813 

4 . 37463E-03 

95.625 

4 . 36005E-03 

98.438 

4 . 07728E-03 

101.250 

3 . 96532E-03 

104.063 

4. 19893E-03 

106.875 

4.33502E-03 

109.688 

4 . 54302E-03 

112.500 

4 . 82269E-03 

115.313 

5 . 15307E-03 

118.125 

5 . 38522E-03 

120.938 

5 . 35189E-03 

123.750 

6 . 11372E-03 

126.563 

6 . 90978E-03 

129.375 

8 . 13378E-03 

132.188 

1 . 15220E-02 

135.000 

1 . 31309E-02 

137.813 

2 . 00602E-02 

140.625 

2 . 26719E-02 

143.438 

2 . 79511E-02 

146.250 

3 . 40878E-02 

149.063 

3 . 29420E-02 

151.875 

2.90356E-02 

154.688 

2 . 12366E-02 

157.500 

1 . 47333E-02 

160.313 

1 . 38721E-02 

163.125 

1 . 02886E-02 

165.938 

2 . 00481E-02 

168.750 

1 . 78127E-02 

171.563 

1 . 67157E-02 

174.375 

2 . 12180E-02 

177.188 

8 . 56301E-02 

180.000 

1 . 84094E-01 

Test  integral  of  phase  function  is:  1.06846E+00 


31 


Literature  Cited 


Abramowitz,  M.,  and  I.  A.  Stegun,  Eds.,  1970,  Handbook  of  Mathematical  Functions.  Applied 
Mathematics  Series  55,  Chaps.  8,  10,  U.S.  Government,  Printing  Office,  Washington, 
D.C. 

Bohren,  Craig  F.,  1983,  Absorption  and  Scattering  of  Light  bv  Small  Particles,  John  Wiley  & 
Sons. 

Lentz,  William  J.,  1976,  Generating  Bessel  Functions  in  Mie  Scattering  Calculations  Using 
Continued  Fractions.  Applied  Optics,  15:3. 

Miller,  August,  1983,  Mie  Code  AGAUS  82.  ASL-CR-83-0 100-3,  U.S.  Army  Atmospheric 
Sciences  Laboratory,  White  Sands  Missile  Range,  New  Mexico. 

Morse,  Philip  M.,  and  K.  U.  Ingard,  1968,  Theoretical  Acoustics.  McGraw-Hill  Book 
Company. 

Nguyen,  X.  Nguyen,  1993,  Conversion  of  AGAUS  Single  Precision  Code  to  Quadruple 
Precision  Code.  Contract  No.  DAAD0791M0664,  U.S.  Army  Research  Laboratory, 
White  Sands  Missile  Range,  New  Mexico. 

Noble,  John  W.,  and  Harry  J.  Auvermann  (Draft),  1993,  The  Effect  of  Large  and  Small  Scale 
Turbulence  on  Sound  Propagation  in  the  Atmosphere.  Battlefield  Environment 
Directorate,  U.S.  Army  Research  Laboratory,  White  Sands  Missile  Range,  New  Mexico 
88002-5501. 

Van  de  Hulst,  H.C.,  1957,  Light  Scattering  ov  Spherical  Particles.  John  Wiley  &  Sons,  Inc. 


33 


Appendix  A.  Notes  on  Acoustic  Scattering  by  a  Spherical  Turbule 

Dr.  G.  H.  Goedecke 
Department  of  Physics 
New  Mexico  State  University 
Las  Cruces,  New  Mexico  88003 

The  problem  of  acoustic  scattering  by  a  spherical  turbule  having  uniform  reference  density 
different  from  that  of  the  surroundings  but  having  no  flow  velocity  is  solved  in  these  notes. 
Complete  analytic  solutions  for  the  differential  and  total  cross-sections  for  an  ideal  fluid 
(no  viscosity)  are  given.  The  analytic  expressions  can  be  evaluated  easily  if  subroutines 
are  available  that  generate  values  of  spherical  Bessel  and  Neumann  functions  (jf(x),  n,(x)) 
for  given  arguments  x  and  that  calculate  values  of  the  Legendre  polynomials  P,(ji)  for  given 
argument  n  for  integers  t  >  0.  These  subroutines  exist  in  any  Mie  scattering  code.  In 
the  following,  bold  symbols  such  as  r  are  three  vectors. 

1.  General  Equations  and  Boundary  Conditions 


Figure  A-l.  Constituents  of  the  propagation  region. 

Referring  to  figure  A-l,  imagine  a  massless  but  perfectly  insulating  membrane  separating 
the  two  regions.  Then  the  pressure  must  be  continuous  across  the  membrane;  otherwise, 
the  massless  membrane  would  experience  infinite  acceleration.  If  the  pressure  is 
continuous,  then  the  massless  membrane  will  not  be  deformed. 

Then,  in  the  reference  condition  with  no  waves  present,  we  have 


P\o  ”  nio*#^io  *  Pio  -  "w^a^a > "  P» 


(A-l) 


35 


*10  ■  "ioVio  *  *20  -  'ho*rr20a  P» 


(A-l) 


where  ni0,  n20  are  the  reference  number  concentrations  in  regions  1  and  2,  respectively;  p10, 
P20,  Po  are  pressures;  T10,  are  absolute  temperatures;  and  kg  is  the  Boltzmann  constant. 
Thus,  if  the  reference  number  densities  are  different,  then  the  temperatures  must  be 
different. 

Now,  consider  small  adiabatic  disturbances  away  from  these  reference  conditions.  Then, 
we  have 


Pi  -  /v*>,To  ;  p2  -  pAI* 20* 


(A-2) 


where  y  is  the  ratio  of  specific  heats.  Now,  let 

"1“  BioO  *  «iOvr»  ;  Hj  -  n^a  ♦  (A-3) 

Then,  equations  (A-2),  (A-3),  and  continuity  of  pressure  across  the  boundary  imply 

€j  -  «j  (A-4) 

at  the  boundary,  or  e(r,t)  is  continuous  across  a  boundary. 

In  each  region,  we  have  the  equations 

—  ♦  V  •  (nr)  *  0  (A-5) 

dt 

^  ♦  *  •  V  »  -  Vp  (A-6) 

dt  Mn 


where  M  =  mass  of  molecule. 
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Equation  (A-5)  implies  that  n,0  v±  is  continuous.  To  first  order  in  small  quantities, 
regarding  v  as  small,  this  requires 

"iovu  ■  *»vu  (a-7) 


at  the  boundary,  or  n^.  is  continuous  across  the  boundary. 

Now,  iinearize  equations  (A-5)  and  (A-6)  using 

p  B  p0Q  *  y  «)  (A-8) 


in  each  region.  Then  equations  (A-5)  and  (A-6)  yield 


de 

dt 


*  v 


» -  0, 


(A-9) 


dv 

dt 


-£°L  Ve  >  - c 2  Ve  , 
Mn0 


ci  _  ypo  _  vkBTo 
Mn0  At 


(A- 10) 


Putting  in  a  time  dependence  e'"1,  we  have 

e  =  J-V*v  (A-l  1) 

(CO 


V 


(A- 12) 


Now,  from  equation  (A-7),  with  h  =  outward  unit  normal  to  boundary  surface  at  any 
point  on  the  surface,  we  have,  with  equation  (A- 10) 
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(A- 13) 


«10vu  -  n  •  V«  -  -M-  A  ■  Ve,, 
10  u  ia  1  iu>M  1 

njjVj  *  •  v«,  =  -^L  if  •  Vej. 

™ 2i  <«  ^  tuAf  ^ 


From  this,  we  see  immediately  that 


if  •  v«j  ■  if  •  V«j 


(A- 14) 


at  the  boundary,  or  &  •  Ae  is  continuous  across  the  boundary. 

Combining  equations  (A- 11)  and  (A- 12),  we  have 

A2e,  *  k2e1  ■  0  ,  A2ej  ♦  mVe 2  *  0,  (A- 15) 


where 


k 


(A- 16) 


and  m  =  "refractive  index"  of  region  2  with  respect  to  region  1. 

Thus,  we  have  to  solve  equation  (A- 15),  subject  to  the  boundary  conditions  (A-4)  and 
(A- 1 4);  that  is. 


e  is  continuous  across  the  boundary, 
A'Ve  is  continuous  across  the  boundary. 


(A- 17) 


These  are  exactly  the  same  conditions  and  equations  that  apply  for  potential  scattering  by 
a  uniform  potential  in  quantum  mechanics;  e  here  corresponds  to  the  wave  function  in 
quantum  mechanics. 
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2.  Application  to  a  Spherical  Region 

We  let  region  2  be  a  sphere  of  radius  a  centered  at  the  origin.  We  assume  a  plane  wave 
of  unit  amplitude  incident  in  the  z-direction: 

m 

«0 W  -  -  £<2«  ♦  l)i%kr)P{c(xQ).  (A  18) 

•  i-o 


The  second  equality  expresses  the  plane  wave  in  terms  of  the  spherical  Bessel  function 
j,(kr)  and  the  Legendre  function  P,(cos0).  So  we  write 


e^r.0)  -  £(2«  ♦  l)i%(kr)  *  atf\kr)}P, (cosB) 

4-0 


(A- 19) 


«j  (r,0)  =  (2t  +  l)itb^l(mkr)Pl(cosQ),  (A-20) 

«-o 


as  the  general  solutions  of  equation  (A' 15)  in  regions  1  and  2,  respectively.  The 
coefficients  (a,  ,b,)  are  to  be  determined  using  the  boundary  conditions  (A- 17),  which  for 
a  sphere  are 


ejfa.e)  =  ej(a,6)  (A-21) 

~(«,0)  =  ~(^0)-  (A-22) 

dr  dr 


The  scattered  wave  is  the  second  term  in  equation  (A-19).  For  r  -*■  oo,  we  know  that 


h?\kr)  -  <-*) 


(A-23) 


so  the  scattered  wave  is 
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(A-24) 


e,(r,Q)  -  —  •  ^-£(21  +  l)affizx&), 
r-*>  r  «,«o 


whereby  the  scattering  amplitude  f(k,0)  is 


Am  -  if  (21  *  lV^coae). 
<*(-0 


The  differential  scattering  cross-section  is 

°(m  -  IM9)I  2- 

The  total  scattering  cross-section  may  be  obtained  in  two  ways. 
First, 

<*,(*)  -  ~E  E  &  *  W*  + 

^  i  ✓ 

But  the  integral  is  equal  to 


2<M 


where  8,r  is  the  Kronecker  delta.  So 


°JW  =  ~E(2f  +  Dl  °i  I  2- 

At  i-o 


The  optical  theorem  yields  the  extinction  cross-section 


(A-25) 


(A-26) 


(A-27) 


(A-28) 


(A-29) 
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(A-30) 


o,(*)  »  ^SSm(A 0'))  -  ^jE(2I  * 

*  It  f-o 

a.(*)  =  -^£<21  *  l)**<g 
t2  w» 


since  P{(cos  0°)  =  1,  for  all  t,  and  3m(-ia,)  =  3m(-i9te(a,)  +  3m(a,))  =  -9?e(a,). 
If  the  refractive  index  of  region  2  is  real,  as  it  is  here,  then  we  should  have 


a.(k)  =  o,(k). 


(A-31) 


So,  once  we  have  found  the  a,,  we  can  evaluate  f(k,0),  cr(k,0),  and  a,.  In  order  to  find  the 
af,  we  apply  the  boundary  conditions  (A-21)  and  (A-22)  to  equations  (A- 19)  and  (A-20). 
This  yields 


jt(ka)  *  aji,m(ka)  =  bjt(mka),  (A-32) 

j,'(ka)  *  afi™\ka)  =  mb/'(mka).  (A-33) 


Here,  the  prime  means  "derivative  with  respect  to  the  argument.  ” 
Solution  of  these  equations  for  a,  yields, 

1  jt(mx)fitw'(x)  -  mj'(mx)hfw(x) 

-  a  - - - - - - , 

ai  j,(mx]j,'(x)  -  mj,Xmx)jt(x) 


where 


x  «  ka. 


(A-35) 


Since  h,(l)  =  j,  +  in,,  we  have 
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I  ,  _{1  +  ~  mj,(mx)n{x)] 

at  j,(mx)jt'(x)  -  mj'(mx)jt(x)  ’ 


-  -  1  -  r, 


where  Tt  is  the  ratio  above,  and  is  real  for  real  m.  Thus 

, _ i_  _  (l  - 

*  i  ♦  «r,  i  ♦  r,1  ‘ 


Therefore, 


l  ♦  r  * 


From  equations  (A-38),  (A-30),  and  (A-29),  we  see  that  indeed 


i-o 


i  +  r/ 


where 


r  B  jfi> «)w/(x)  -  iwy/(mx)n,(x) 
1  ;,(mxV/W  -  rnjt'(mx)jt(x) 


Depending  on  m  and  x  =  ka,  T,  -*  0  for  some  £  >  that  is,  we  can  neglect  £ 
The  larger  m  and  the  larger  ka,  the  larger  Roughly,  *  mka. 

computational  purposes,  we  take  all  sums  from  £  =  0  to 

For  example,  suppose  mka  <  1.  Then  the  dominant  term  is  £  =0,  and  we  get 


(A-36) 


(A-37) 


(A-38) 


(A-39) 


(A-40) 


So  for 
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.  o(*,0) 


1 


(A-41) 


Mfi) 


j_  . _ i_ 

a  i  +  «r0 


*?(  i  *  r0J) 


independent  of  0;  so  the  scattering  is  especially  symmetric  (monopole).  This  is  to  be 
contrasted  with  the  lowest  order  multipole  for  the  corresponding  electromagnetic  wave 
problem,  which  is  dipole  in  character,  with  a(k,0)  dependent  on  B. 

For  this  case,  mka  <  1,  we  can  evaluate  T0  easily.  We  need  only  the  asymptotic  forms 
of  (jo.  °o)  for  small  argument. 

From  equations  (A-12)  and  (A-24),  we  see  that  the  scattered  wave  velocity  field  is  given 
in  the  radiation  zone  (r  -*  «)  by 


— vc,(r,e>  —  —  f  cflm. 

!(■)  f--4«  T 


(A-42) 


Thus,  the  fluid  velocity  of  the  scattered  wave  is  radial  (longitudinal).  This  is  to  be 
contrasted  with  the  electric  and  magnetic  fields  in  the  scattered  wave  in  an  electromagnetic 
scattering  problem,  which  are  transverse  to  f . 

At  any  rate,  we  have  a  complete  analytic  solution  to  the  problem  of  scattering  of  acoustic 
waves  by  a  sphere  of  uniform  reference  density  different  from  that  of  the  surroundings. 
Merely  evaluate  the  T,,  £  =  0,  using  standard  algorithms  for  the  spherical  Bessel  and 
Neumann  functions  (j,,  n,);  then  substitute  the  results  into  equation  (A-37)  for  a,,  then 
substitute  a,  into  equation  (A-25)  for  f(k,0)  using  standard  algorithms  to  evaluate  the 
Legendre  functions  P,(cos0). 
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Appendix  B.  The  Four  Size  Distribution  Models  of  ASCT 


Type  0:  Arbitrary  User-Supplied  Discrete  Model 

This  distribution  type  allows  the  user  to  enter  from  1  to  513  pairs  of  data.  NRAD1 
is  the  pairs  of  radii  and  absolute  number  densities.  The  radius  must  be  entered  in 
meters  and  the  absolute  number  density  in  particles  per  cubic  meter  per  meter.  The 
input  file  has  the  following  format:' 

NWAVE  NIDSTP  ICQ  NANG  IANG  NBINS  IEO  NEOU  NUN1T  MQRTE 
DELTA  ATOP  (if  IDELT  =  1) 

IDSTP  NRADI  0.  0.  0.  0.  0. 

R1  F(R1) 

R2  F(R2) 


RNRADI  F(RNRADI) 

WAVE  EM  CAY  0.  0.  0.  0.  0. 

Type  1:  Log-Normal  Distribution  Model 

This  model  is  used  to  compute  the  cross-sections  and  phase  function  for  a  log-normal 
distribution  of  turbule  sizes  and  is  described  by  the  following  formula: 


Ar)  -  =1. 

r  2  ln(o) 


(B-l) 


The  input  file  has  the  following  format: 

NWAVE  NIDSTP  ICQ  NANG  IANG  NBINS  IEO  NEOU  NUNIT  MQRTE 
DELTA  ATOP  (if  IDELT  =  1) 

IDSTP  RLO  RHI  RBAR  SIGMA  DENS 
WAVE  EM  CAY  0.  0.  0.  0.  0. 

RLO  -  The  minimum  radius  (m)  at  which  the  model  is  to  be 
cutoff. 

RHI  -  The  maximum  radius  (m)  at  which  the  model  is  to  be 
cutoff. 

RBAR  -  The  radius  (m)  at  which  the  relative  number  density 
peaks. 

SIGMA  -  The  standard  deviation  to  be  used  (not  LN(SIGMA)). 

DENS  -  The  number  of  particles  per  cubic  meter. 
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Type  2:  Single  Turbule  Model:  Efficiency  Factors/Cross-Sections 
Phase  Function  Calculations 

Used  to  compute  efficiency  factors  (ICQ  =  0)  or  cross-sections  (ICQ  =  1)  and  phase 
function  for  a  single  turbule  radius.  The  input  file  has  the  following  format: 

NWAVE  NIDSTP  ICQ  NANG  IANG  NBINS  IEO  NEOU  NUNTT  MQRTE 
DELTA  ATOP  (if  IDELT  =  1) 

IDSTP  RADIUS  0.  0.  0.  0.  0. 

WAVE  EM  CAY  0.  0.  0.  0. 

Type  3:  Special  Table  Generation  Mode 

This  is  not  a  distribution  model  but  is  an  operating  model  used  for  calculating  and 
printing  tables  of  the  efficiency  factors  as  a  function  of  radius.  The  input  file  has  the 
following  format: 

NWAVE  NIDSTP  ICQ  NANG  IANG  NBINS  IEO  NEOU  NUNIT  MQRTE 
DELTA  ATOP  (if  IDELT  =  1) 

IDSTP  RLO  RHI  DELR  0.  0.  0. 

WAVE  EM  CAY  0.  0.  0.  0.  0. 

RLO  -  The  minimum  turbule  radius  (m)  to  be  used. 

RHI  -  The  maximum  turbule  radius  (m)  to  be  used. 

DELR  -  The  increment  in  radius  (m)  to  be  used. 
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Appendix  C.  AGAUS  Formulas 

The  formulas  implemented  in  program  AGAUS  for  calculating  the  efficiency  factors  and 
the  phase  function  are  given  below.  They  are  included  here  to  compare  the  corresponding 
acoustic  formulas  for  efficiency  factors  and  phase  function  more  easily. 

Expansion  Coefficients 


A.(mz)  t 

[-Z—1  *  -]8k[5tto]  -  «*K,.,(x)j 
m _ x _ _ 

y4.(mx)  | 

*  -K/x)  -  $,_,(*) 
m  x 


(C-l) 


[mAt(mx)  *  -]**[?,(*)]  - 

b,- - - - ,  (C-2) 

[mAfmx)  *  -  $,_,(*) 


A.(mx)  «=  — — 
mx 


J,.  i(”**) 

y,(«x) 


(C-3) 


Efficiency  Factors 


<?«  =  4E(2?  +  +  M 

Jr  1*1 


(C-4) 


<?„  ■  *  '<1  -.1 !*  I  M  1  (C‘5) 

Jr  i*i 

1  «*  2 

=  4l  E(2«  ♦  1K-DH  +  *,)l  (C-6) 

XT  (.] 
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Scattering  Amplitudes 


St(0)  -  Vi(0)1 


s#>  -  l  W®>+ 

1-1  “*  +  AJ 

Intensity  Factors 

iv(6)  -  I  Sj(0)|  2 

4(6)  -  I  S3(0)|  2 

Recurrence  formulas  for  Pi  and  Tau  Functions 
«,(0)  -  cos0^f^n,.1(0)  -  ”-n,.2(e) 


T|(0)  »  cos0[nt(0)-  it,.j(0)]  -  (2t-l)sin20n<_1(0)  +  jt^O) 


Phase  Function 


_  [<‘,(6)  ♦  4(0)1 

2i?C„ 


For  all  the  above  f  =  1,2,3,.... 


(C-7) 

(C-8) 

(C-9) 

(C-10) 

(C-ll) 

(C-12) 

(C-33) 
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